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We use the experimental technique of thermal modulation reflectometry to study 
the relatively small temperature dependence of the optical conductivity of supercon- 
ductors. Due to a large cancellation of systematic errors, this technique is shown 
to a be very sensitive probe of small changes in reflectivity. We analyze thermal 
modulation reflection spectra of single crystals and epitaxially grown thin films of 
YBa2Cu307^,5 and obtain the atr 2 F(iv) function in the normal state, as well as the 
superconductivity induced changes in reflectivity. We present detailed model calcu- 
lations, based on the Eliashberg-Migdal extension of the BCS model, which show 
good qualitative and quantitative agreement with the experimental spectra. 
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I. INTRODUCTION 



Optical spectroscopy in the infrared region provides a valuable tool to study low-energy 
excitations in high T c superconductors. In principle such measurements can provide infor- 
mation about the superconducting energy gap (if any) and the pairing mechanism. A lot 
of discussion is still focussed on the problem of how to distinguish electronic contributions 
from those which are due to phonons, and to what extent these two are coupled. The latter 
deserves special attention, because in the conventional superconductors electron-phonon cou- 
pling is believed to be responsible for superconductivity. While the role of electron-phonon 
coupling is not very clear in the new high T c materials, there exists a tendency in the theo- 
retical community to consider conventional coupling mechanisms such as electron-phonon as 
a necessary ingredient to obtain pairing. It is not yet clear however whether an exceptionally 
large electron-phonon coupling is really required, or whether a small germ suffices, which 
is then further 'boosted' by other mechanisms such as interlayer pairing (as was recently 
proposed by P.W. Anderson [|l]). This uncertainty emphasizes the importance of a proper 
understanding of electron-phonon coupling in these materials if we want to make progress 
towards a theory of high T c which is of practical importance, e.g. able to predict T c for new 
compounds. Although much data of high quality has been obtained on a wide variety of 
cuprate high T c superconductors, it has turned out difficult to obtain a clear interpretation 
of these data. In particular there have been many papers where the observation of a 'clean' 
gap has been claimed, but closer inspection reveals, that there is a considerable amount of 
residual absorptivity at energies below the presumed gap. Some of the difficulties are con- 
nected with experimental problems concerning the complicated structure of such systems, 
but others are related to the theoretical interpretation of the results obtained. Two exam- 
ples are that the value these gaps would have is far too high to be compatible with a simple 
BCS picture, and the position of these features doesn't appear to shift to zero frequency if 
T — > T c . Instead a gradual filling of the gap region is observed with a change of slope of the 
intensity in this region at the phase transiton, which can be fitted in a phenomenological 
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way using a two-fluid description of the superconducting state. 

In the conventional theory superconductivity is a consequence of electron pairing me- 
diated by the exchange of the Bose-like excitations. This pairing leads to the appearance 
of an energy gap and hence to zero absorption for frequencies less than 2 A. In the optical 
data it is possible to determine this threshold by measuring the departure of the reflectivity 
from unity, or of the absorptivity from zero. The most informative value is the optical con- 
ductivity (a). But in order to calculate a it is necessary to determine with great accuracy 
small departures from 100% reflectivity, followed by a Kramers-Kronig analysis to obtain the 
phase of the reflectivity. When the absolute reflectivity is close to 100% the progression of 
experimental errors in this procedure often leads to large uncertainties in the conductivity. 

For experimental reasons 2A was determined in conventional superconductors as a max- 
imum of a ratio of the reflectivity in the SC state to the normal state. For HTS there are 
several difficulties when one tries to apply this procedure due the high critical temperature 
T c itself. Between zero temperature and T c some intrinsic changes in the electronic subsys- 
tem occur, due to which we have practically two different substances at low (helium) and 
normal state temperatures. This is more important if the mechanism of pairing is of an 
electronic nature. 

II. THERMAL MODULATION REFLECTOMETRY 

Another method of the determination of the energy gap has been used by Abel et al. @ . 
It provides a possible way to avoid some of the afore mentioned difficulties (e.g. problems 
with the absolute value of reflectivity and Kramers-Kronig transformations). Abel et al. de- 
termined the ratio of optical reflectivities at two close temperatures R(u, T)/R(u, T + ST), 
and the maximum of this ratio was ascribed to the superconducting energy gap 2A(T). 
This feature reflects the fact that the main change below T c is due to the decrease of A in 
the single-particle excitation spectrum. However, in real superconductors the situation is 
more complicated. First, for T > the gap itself is not well defined. Second, we have a 
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strong temperature dependence of the number of thermally excited Bose-like quasi-particles. 
As we shall see below, this is the dominating effect, and gives us a possibility to extract 
information about the spectrum of these intermediate bosons. The third complication is 
connected to the so-called Holstein shift. The optical conductivity contains contributions 
from the intermediate bosons which are responsible for the pairing because they are coupled 
to the electron-hole excitations (boson- assisted conductivity, see e.g. HU). The difference 
between the normal and the superconducting state is the following: A threshold point of 
an absorption which in the former case takes place at the characteristic boson frequency 
Qo should in the latter case be shifted to Qq + 2 A. This effect was observed in conven- 
tional superconductors and played an important role in establishing the phonon-mediated 
nature of superconductivity in these systems However, the absence of such a shift in 

the reflectivity spectra of high T c superconductors is one of the most important objections 
against the conventional mechanism of superconductivity in these compounds J7|||. The 
thermo-modulation method, in which the ratio of reflectivities at two close temperatures is 
analyzed |2|], gives the possibility to observe the boson anomalies in the normal as well as 
in the superconducting state. In the present paper we make a detailed comparison between 
model calculations based on strong-coupling theory, and experimental thermo-modulation 
reflection spectra. The data used in the present analysise are in good agreement with those 
obtained by Abel et ai, but cover a wider range of frequencies and temperatures. 
The main advantage of measurements of R(u,T)/R(u,T + ST) between two close temper- 
atures is that all thermal effects connected with extrinsic factors (e.g. the experimental 
setup, spurious signals) are compensated and the ratio reflects only the temperature depen- 
dence of the occupation numbers of the bosonic and fermionic excitations. If the change of 
temperature 5T <C T, it is possible to expand 



From the expression for the reflectivity at normal incidence R = |(1 — y / e)/(l + v^e)] 2 one 
obtains the following exact expression for the thermo-reflectance coefficient 




(1) 
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r(u,T) = -din R(u,T)/dT = -2Re 



de/dT 



(2) 



Although in the comparison which we will make to experimental optical data we will use 
numerical calculations based on the full solution of the Eliashberg equations, we will make 
some further approximations in order to demonstrate that r(u, T) reflects the af r F function. 
We first notice, that for a good metal in the frequency region under consideration (u <C uj p ) 
\e'\ ^> 1, so that we may write r(u,T) = ARed(e~ 1 ^ 2 )/dT. As can be seen from Fig. 8 of 
Ref. 0.25 < e"/e' < 0.4 for the relevant frequency range in the superconducting state. 
Although in the next section we will make a comparison between experiment and theory 
using the exact expression for the thermo-reflection (Eq. 2), in the present qualitative 
discussion we only consider the 'clean' limit, where (e"/e') = {ujt)~ 1 may be treated as a 
small parameter in a Taylor series expansion. The leading order of this expansion is 

d e" 



r(u } T) 



(3) 



dT | e' | 3 / 2 

For the description of the dielectric function in the normal state we use the extended Drude 
expression 



e(u,T) 



(4) 



m*(uj)/m + i/iuriuj)) 

We will make the approximation here, that the main temperature dependence enters through 
the parameter t(oj). As we consider here the range of frequencies where lot ^> 1, we obtain 
the simple expression 

2 dr- 1 



r(w,T) 



oo p dT 



(5) 



As is shown in the Appendix, for sufficiently low frequencies the temperature derivative of 
the optical scattering rate is proportional to the transport spectral function, whereas at high 
frequencies it is proportional to X tr , so that 

Xtr (T > fi ) 



r(u,T) 



47T 



2n 2 T „,2 
3u> 



al(u)F(u) (T « fi ) 



(6) 



So we see, that thermal modulation reflectometry can provide rather direct experimental 
information on the transport spectral function. 

In the superconducting state it is possible to obtain a useful expression in the region where 
— e' ~ ^f^r ^> 47rCT1 ^ aJ ' T ) (i.e., when R ~ 1), where o\ = cue" /(An) is the real part of the 
optical conductivity. Here Al is the London penetration depth. In this case 

r(u,,T)=Snc-^W- d -^l (7) 

In the Appendix we show that in the superconducting state the dominant term in 9ai ^' T ' > is 
proportional to a^ r F(u — 2A). As a result we obtain, disregarding the slowly varying term 

r(u, T) = ^f hT aUu - 2A)F(co - 2A) (8) 

At intermediate temperatures the dielectric function can sometimes be approximated with 
the two-fluid interpolation formula e(u,T) = e s f s (T) + e n (l — f s (T)), as was supported 
experimentally for Yl^CusOy.^ and theoretically in the case of strong or intermediate 
coupling ]nj, where e s (co) and e n (uj) are the dielectric functions at T = 0, and at T > T c 
respectively. The function f s (T) is proportional to the number of superconducting electrons 
and is assumed to be of the form f s (T) = 1 — (^) where v is some exponent (usually 
v = 4). For T < T c we see that r(u, T) oc Re[a n (uj) — c s (uj)} and has its maximal amplitude 
directly below T c . 



III. COMPARISON OF EXPERIMENTAL DATA WITH STRONG-COUPLING 

CALCULATIONS 

In Fig. la the experimental values of the thermo reflectance coefficient r(u,T) are dis- 
played for an ab-oriented thin film of YBa2CusOjs with T c = 90K JT^j for temperatures 
between 30K and 140K. The ratios were calculated from data differing by 10 K. We see a 
feature in the frequency region uj < 1000cm -1 both in the normal state and in the super- 
conducting state. We made a similar analysis of the data obtained by Bauer on a single 



7 



crystal of YBa2Cu3C>6.8 with a T c of 81-86 K ||14|| . As only reflectivity data at 10, 60 100, 



and 150K were taken, the ratio's Ta ^ Ti ^[^jl^pl) (displayed in Fig. lb) are calculated for 
the pairs (T 1 ,T 2 ) = (10,60), (60, 100) and (100, 150). Although one expects deviations from 
a pure thermal derivative in this case, the thermal modulation spectra obtained from both 
sets of data are actually quite similar. The origin of the larger value of r between 100 K 
and 150 K in Bauer's data is not clear, and may be due to the lower oxygen-concentration 
with a T c of 80-85 K in this sample. In the normal state all features end at the frequencies 
near 700cm -1 , which corresponds to the range of phonon frequencies for these materials. 
According to expression 6 the shape of r(u),T) should be proportional to the spectrum of 
intermediate phonons a^ r F(u>). This gives us the possibility to speak about the intermediate 
boson contribution to the optical properties of high T c systems. 

Let us first calculate the thermo-modulation spectra using the weak coupling BCS model 
in the clean limit, either assuming a single gap at 220cm _1 (Fig. 2a) or a distribution of 
gaps between and 500cm _1 (Fig. 2b). We see, that the calculated r(u,T) has only one 
feature corresponding to the energy gap. Also the calculated thermomodulation effect on 
the reflectivity is much larger than the experimental values. 

To obtain a better understanding of these features we carried out model calculations 
of the reflectivity assuming an o?F function with a single broad peak at ujq = 350cm _1 , 
the constant of interaction being A = 1.5, and the bare uj p = 3eV, which gives a critical 
temperature T c = 87K. This shape of a 2 F{oj) was used in |H| These parameters lead to 
a linear dependence of the resistivity in the normal state, with a slope that corresponds to 
experimental values fl5| . To match the 164K data it is either necessary to take or a larger 



ujq or larger A (~ 2.5). The results are shown in Fig. 3. 

According to expression 8, for T ^ T c the spectrum of the intermediate bosons should 
be shifted by 2 A. The experimental data indeed show such a shift and make it possible to 
estimate the value of 2A (if A& is actually a distribution of gap- values due to e.g. anisotropic 
pairing, the shift corresponds to a gap-value averaged over the Fermi-surface) to be 250 to 
300 (cm -1 ), corresponding to a ratio 2A/T C w 4 to 5. It is interesting to note the negative 
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contribution to r(u,T) above 2A + Q, which is a consequence of the modification of the 
optical conductivity (second term in Eqffi]). At intermediate temperatures the ratio r(u,T) 
behaves according to the two-fluid model and reaches a maximal amplitude just below T c . 



IV. CONCLUSIONS 

We demonstrate that thermal modulation reflectometry can be used to record small 
changes in reflectivity of superconductors, even if the reflectivity itself is close to 100%. We 
observe that the thermal modulations of superconducting YBa2Cu 3 7 _s are quite small, but 
well reproducible from sample to sample. Although the observed effects are much smaller 
than a weak coupling BCS-type calculation in the clean limit, qualitatively good agreement 
with a strong-coupling type model calculation is obtained, where A = 1.5, and u p = 3eV is 
assumed, and an a 2 F function with a single broad peak at 350 (cm -1 ) is taken. 
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VI. APPENDIX 



A. Normal state 



Let us first consider the reflectance in the normal state. The temperature dependence of 
the optical scattering rate can be expressed using the relation of Ref. |10[ (valid for lot ^> 1) 



(9) 

I Jo 00 dn a l(n)F(n) [200 coth (§) + (w-n) coth (*=pj -(u + n) coth (ago) 

where a^ r (u})F(u) is a transport spectral function of the electron-boson interaction 

9 



al{u)F{u) 



N(0) 
4v F 



<| M kk , | 2 (V k - V kl f8{n % _ %l - u) » 



(10) 



where n and M are an intermediate boson frequency and a matrix element, and <C .... ^> 
denotes an average over the Fermi surface. In the limiting cases, where the temperatures 
are either much smaller or much larger than the characteristic frequency Qq, we have 



t(u,T) 



-i 



4ttT / °° d(\nn)al(VL)F(n) = 2n\ tr T 



2 - r dnal(n)F(n)(uj - n) + f£T 2 aiMFM (t « no) ^ 

The first term has been previously obtained by Allen M. The first derivative with respect 
to temperature becomes 

2n\ tr (T > n ) 



8t{uj,T)- 1 
df 



47r 3 T ^,2 
3lu 



(T«n c 



(12) 



B. Superconducting state 



For T < T c the expression for the boson-assisted conductivity was obtained by Allen 



at(uj,T) 



6u, 2 



<| M kk , | 2 (\4 - ^fc')Vfc(l - /fc') (l 



E k E k> 



5{E k + + ntt, - u) > 



(13) 



Here = y ef. + A| is a quasiparticle spectrum, is the superconducting gap and f k is the 
Fermi-distribution. For low temperatures T A it is possible to neglect the temperature 
dependence of the gap, and we obtain 



dai(uJ,T) __ / vre 2 2ir 2 nT 
dT 1 6u> 2 m 



(14) 



{-AaiF{^ - 2A) + ;;- 2A dnaiF(ny A{u _ Q _ A) 

The first term is proportional to the Eliashberg function af r (uj)F(uj) shifted by 2A (Hol- 
stein shift). The second term is weakly varying with frequency and is a consequence of a 
modification of the conductivity at uj > 2A + Q. 
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FIGURES 

FIG. 1. a: Experimental r(u,T) with E _L c for an epitaxial thin film of YBa2Cu307_,5. From 
top to bottom: T= 35, 55, 75, 95, 115 and 135 K. The curves have been shifted vertically with 0, 
-0.0005, -0.0010 FT 1 etc. 

b: The same for a single crystal. Temperatures are 35, 80 and 125 K. Vertical offsets have been 
given of 0, -0.001 and -0.002 K _1 . 

FIG. 2. a: Theoretical calculation of r(u,T) based on the BCS model (2A = 3.5ksT c ) using 



the method of Ref. 16]. Temperatures are 35, 55 and 75 K from top to bottom. Vertical offsets 
have been given of 0, -0.001, and -0.002 K _1 . 

b: Theoretical calculation of r(u, T) based on a BCS-like model with a gap distributed between 
and 8ksT c . Temperatures and offsets are as in Fig. 2a. 

FIG. 3. Theoretical calculation of r(u,T) based on the strong coupling formalism with param- 
eters as explained in the text. From top to bottom: T= 35, 55, 75, 95, 115 and 125 K. The curves 
have been shifted vertically with 0, -0.0005, -0.0010 K" 1 etc. 
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